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1— h Abstract 

FAST (Fusion Advanced Studies Torus), the Italian proposal of a Satellite Facility to ITER, is a 

CN 

compact tokamak (Rq = 1.82 m, a = 0.64 m, triangularity 5 = 0.4) able to investigate non-linear 

a 

dynamics effects of a-particle behavior in burning plasmas and to test technical solutions for 

O 

the first wall/divertor directly relevant for ITER and DEMO. Currently, ENEA is investigating 
r j^ 1 the feasibility of a superconducting solution for the magnet system. This paper focuses on the 

analysis of the CS (Central Solenoid) magnet thermal behavior. In particular, considering a 

superconducting solution for the CS which uses the room available in the resistive design and 
c/3 referring to one of the most severe scenario envisaged for FAST, the heat load of the CS winding 

c/3 pack due to AC losses is preliminarily evaluated. The results provide a tentative baseline for the 

definition of the strand requirements and conductor design, that can be accepted in order to fulfil 

the design requirements. 
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I. INTRODUCTION 



FAST (Fusion Advanced Studies Torus) is the Italian proposal for a new European satel- 
lite tokamak reactor aimed at supporting ITER activities and anticipating some DEMO 
relevant physics and technology issues [HE]- It has been conceived as a compact (Ho = 1.82 
m) machine working at high field (BT up to 8.5 T) and high plasma current (Ip up to 8 MA). 
Currently, FAST magnetic system is designed with 18 TF (Toroidal Field), 6 CS (Central 
Solenoid) and 6 PF (Poloidal Field) resistive coils, cooled by helium gas flow at 30 K. A 
feasibility study to verify whether a superconducting (SC) solution for the whole magnetic 
system would be possible or not, avoiding any major modifications to the machine geometry 
or scenarios, is currently under study at ENEA |3J. This paper focuses on the proposed 
superconducting solution for the CS magnet in order to assess its feasibility. Since high 
operative currents together with both high fields and high field rates are envisaged during 
the operation of the CS, the evaluation of the AC losses heat loads and the corresponding 
thermal behavior of the SC cables are mandatory. Therefore, starting from the layout pro- 
posed in [3J, based on the 7.5 T H-mode reference scenario of FAST [fl [2], the heat load of 
the CS superconducting cable originated by AC losses is evaluated in order to provide both 
a tentative baseline for the definition of the strand and cable requirements and the inputs 
needed for a future, complete thermohydraulic study. 

II. DEFINITION OF THE PROBLEM 
A. Design details 

The SC feasibility study, presented in [3J and used in the present study, has been focused 
on the H-mode reference scenario (BT = 7.5T, Ip = 6.5 MA, other engineering parameters 
available in pQ or [3]), which is one of the most challenging scenarios among those foreseen 
during FAST activities. It is worth noting that this scenario is characterized by a relatively 
short duration, having a flat top lasting for ~ 12 s only. The time evolutions of the operative 
currents I op for the CS1 modules only (see later) are reported in Fig. [lj 

As it can be seen from Fig. [TJ the H-mode reference scenario can be divided in the 
following main phases: the charge of the coil, lasting 8 s; a fast ramp down of the field B 
(with a very large B in the initial 0.04 s), lasting ~8 s; the flat top, lasting ~ 12 s; the final 
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FIG. 1: Operative current vs time, in CS1 modules 



discharge of the coil, lasting other 8 s. 

Based on the room available in the resistive design, the SC version of the CS coil has been 
conceived. It consists of a rectangular winding pack (WP), inserted into a stainless steel 
casing, wound from a rectangular conductor with an aspect ratio lower than 2, having long 
twist pitch values, low void fraction (~ 29%) and no central channel. These choices, mainly 
aimed at giving better load support to the NbsSn strands and consequent limited perfor- 
mance degradation, derive by many measurement campaigns and studies performed in the 
last years [5]. Note that the choice of M^n for the SC strand is required by the maximum 



magnetic field (around 17.7 T, see Section II B) arising in the coil. The details of the design 



for all the 6 modules making the CS coil are reported in [3|; since the central modules (CS1U 
and CS1L) experience the more stringent conditions (highest field intensities), we limit our 
analysis to them. Their main characteristics, together with details of the WP, are reported 
in Table [H 



B. Magnetic fields 

Another fundamental ingredient of the present analysis is the magnetic field intensity 
experienced by the SC cable. A complete 3-D electromagnetic analysis of the whole scenario, 
including all the field contributions, has been performed with the Opera-3D code. All the 
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Cable inner dimensions (mm 2 ) 


28 x 35.2 


Cable jacket thickness (mm) 


2.5 


Corner radius (mm) 


250 


Strand diam. (mm) 


0.81, [13] 


A=(l+Cu:nonCu ratio) -1 


0.5 


Cabling pattern 


3x3x5x5x6 


# strands (all SC) 


1350 


VF 


about 29 % 


lo P (kA) 


40.2 


l^peak (T) 


17.7 



TABLE I: central modules CS (CS1U and CS1L) and WP main features [3]. 

three magnetic field components have been computed in function of the z coordinate (the 
CS axis is chosen as z axis, with the origin in the equatorial plane separating the CS1L and 
CS1U modules) and of the radial coordinate (with full computations for inner, outer and 
middle radii of the CS module; in-between values has been linearly interpolated). The time 
evolution has been described by linearly interpolating the values fully computed at selected 
instants of the transport current rate changes. The maximum peak field 17.7 T is found at 
z=+12 cm from the equatorial plane. In the following we consider the pancake (belonging 
to the CS1U module) subjected to this peak field: being near the equatorial plane, the 
main contribution to the magnetic induction B arises from the z-component, whereas the 
azimuthal one is zero and the radial one is small and hereafter neglected. Therefore, the 
whole pancake cable is assumed to be subjected to a unidirectional transversal field B z . 
Each turn of the pancake has been considered to be subjected to a uniform field equal to 
the one actually present in its axis, i.e. evaluated at its median radius. 

C. Heat loads - AC losses 

The main heat load on the superconducting cable is originated by the applied time- varying 
electromagnetic fields, which cause power dissipation within the superconductor material as 
well as in the normal metal matrix and cable jacket. Neutron heating is neglected in CS 
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coils since it is essentially fully screened by the straight legs of the TF coils surrounding the 
plasma chamber. It is well known [B] that dissipation in SC cables submitted to time- varying 
fields and transport currents consists of three main contributions: hysteresis losses, coupling 
losses, and the so-called "dynamic resistance" originated losses. 

The hysteresis losses Qh occur in Type-II superconductors because of the irreversible 
magnetization, due to the pinning of vortices, which arises when the external magnetic field 
B is swept. The main parameter governing these losses is the penetration field B p here 
defined as the field at which an applied increasing field fully penetrates the cable, starting 
from the virgin condition. The penetration field B p depends on the geometry, dimensions 
and critical current densities. It is also worth noting that B p depends on the temperature 
T and on the external magnetic field through the critical current density. 

The coupling losses Q c are due to coupling currents which arise in multifilamentary wires 
and multistage cables, flowing in the superconducting wires and closing through the nor- 
mal metal (the matrix and the stabilization layer), where Joule dissipation actually occurs. 
These current are excited by the time-varying flux of the magnetic field threaded to the 
multifilamentary wire. The main parameter governing coupling losses is a characteristic 
coupling time constant r, which dictates the time scale of coupling current transients, and 
which depends on the strands/wires twist pitch and on the inter-filament /strand transverse 
resistivity. Actually, multistage cables introduce multiple (and corresponding volume 
fractions) for each stage of the twisting; here, as often done for simplicity sake, we consider 
a single effective time constant. 

Finally, the "dynamic resistance" losses Qd occur whenever the current carrying cable 
is subjected to a time-varying field: the latter would change the transport current which 
has to be maintained by the current source. The work, and thus the dissipated energy, is 
then provided by the current source itself instead from the external magnetic field, as it 
happens for the two losses above discussed. This latter contribution is often incorporated 
in the hysteresis loss term, since they share a common nature. 

In addition, one has to take into account losses due to eddy currents Q e in pure bulk 
normal metals: in our case, since the strands are all SC, the only normal metal is given by 
the cable jacket. 

According to [0] (and to [7] regarding the eddy current losses), the expressions for the 
power dissipation per unit volume (W/m 3 ) in the static regime (slow field changes, see later) 
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are: 



P c = nr 



(1) 



Ph 



B e 



(2) 



Pd 




(3) 



P e 




2 



(4) 



where the moduli ensure that the power losses are always positive, as they should be, for 
both increasing and decreasing fields; % = It/I c is the transport current normalized over 
the critical current, [IB], n is a geometrical factor (in general taken as equal to 2, more 
often incorporated in tit as a single global quantity), p jacket and l\ are the resistivity of 
the jacket and its outer dimension, normal to the field, respectively. The quantities Bi and 
B e = H/fio are the induction field inside the cable and externally applied field evaluated 
inside the cable, respectively; the dot denotes the time derivative. In the static regime, the 
external and internal fields are related through the following expression [8]: 



The expression for hysteresis losses assume that the actual external field swing is ^> B p . 

As far as geometry is concerned, the system can be described in terms of (isolated) 
cylindrical strands in a transversal applied field: the corresponding penetration depth is B p = 
(l/ir)fi J c d e ff, where d e ff is the effective diameter of the superconducting filaments within 
the strand, and J c is the critical current density of the superconductor (I c = J c A nonCu = 
J C XA, where A and A non c u = AA are the total and nonCu strand sections, respectively). 
The numerical prefactors in Eq. (|2| and ^ have been taken within the same geometry 
(isolated cylindrical strands in a transversal field). For a line of cylindrical cables/strands 
close-packed along the field direction, the slab geometry (with the corresponding geometrical 
factors) could be a possible description [6 J . Since the strands considered in the following have 
a significant thickness of the stabilizing matrix (« 0.2 mm), the filamentary zone within the 
individual SC strands (thickness ~ 0.4 mm) can be considered sufficiently separated (~ 0.4 
mm) in order to be satisfactorily approximated by the isolated cylinder geometry. 

Finally, in order to obtain the power loss per unit length (W/m), Ph and Pd must be 
multiplied by A nonCu = AA, P c too (provided that nr is correspondingly defined as to be 



B e — Bi = T Bi 



(5) 
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referred to the A non c u only), and P e by an effective jacket cross section Aj ac k e t = (h * h — 

11* hi 'id m- 

As initially stated, the above model holds in the static regime, i.e. slowly varying fields 
with respect to the system characteristic time. With rapidly time-varying fields, the screen- 
ing currents in the outer layers of the cable reach J c ; by further increasing B e , the thickness 
of the screening layer increases until all the superconducting filaments reach J c . No further 
screening is then possible and the cable behaves like a solid monolithic conductor. A satu- 
ration parameter is introduced jU] as f3 = tB^/ B PjStrand , where B p>strand = (l/n)fi D e ff\J c is 
the penetration field of an "equivalent" strand, i.e. of a superconducting cylinder of diame- 
ter D e ff having the same J c and magnetization of the cable. As already noted for B p , it is 
worth noting that both B PtStrand and (3 depend on Bi and T through J c . The static regime 
occurs at /3 <C 1, whereas full saturation occurs at /3 ^> 1. In the latter case, expressions for 
power losses are [91 [10] : 

4 2 l-^el 

Pc = ~^~B P strand {l — i ) (6) 
P h = (7) 

P — — R i2 \Be\ fo\ 

1 d J - J v, strand 1, \° ) 

where, in this limit, the following relation holds: 

B e -B % = B p>strand (l - i)sign(B e ) (9) 

For completeness, we mention that in the intermediate regime (given approximately by 
/3 < 1 — i), models are definitely cumbersome: only the model of Ref. [H] allows to explicitly 
calculate all the quantities without complicated numerical computations, but on the other 
hand it fails to describe accurately the crossover between the two limiting regimes [TT] . 



D. Cooling system 

The cooling path is taken in its simplest form: a double pancake cooling is considered 
(although, given the small length w 42 m of a single pancake, also an hexapancake assembly 
could be devised). Therefore, the study can be focused on a single pancake (half of double 
pancake). Typical values for the helium flow are considered: inlet temperature 4.4 K, 
pressure 6 bar and helium flow 5 g/s. 



7 



III. SIMULATION RESULTS 



We compute the AC losses in the pancake located where the applied field B e is maximum. 
In order to simplify the computations, coherently with the scope of the present preliminary 
computations, we take Bi « B e = B (and therefore in Bi « B e = B). Since the considered 
B e (t) is piece-wise linear, in this way a small error is done, involving only small time intervals 
w t after the instants where B e changes. In the static regime, the main parameters of the 
model are the coupling time nr (which governs the coupling losses) and the filament effective 
diameter d e ff (which governs the hysteresis losses). We take nr = 55 ms (related to the 
non-Cu section of the strand), as experimentally determined in [T2] for the orientation of 
the cable with the long side parallel to the applied field, and d e // = 50 jttm as a value 
representative of the filament size in the SC strand here considered. 

The daring choice of all these parameters deserves some further comments: in order to 
meet the very demanding requirements of a large-size CICC (Cable In-Conduit Conductor) 
operating at 17.7 T and 40 kA, the characteristic properties of the most performing Nh^Sn 
strands (non-Cu J c (12 T, 4.2 K, -0.46%)~1500 A/mm 2 [13]) available on the market have 
been considered at this stage of the design, together with conductor design parameters that 
have in principle been optimized for TF magnets operating conditions. Very well conscious 
of the limiting aspects of the present choice for a CS magnet, the study presented here 
should provide a baseline for the definition of the strand requirements (maximum d e ff) 
and conductor design (maximum cable nr), that can be accepted in order to fulfil the 
present design requirements. The final, optimized choices will need to be made as a trade- 
off between minimizing the heat load due to AC losses, while maintaining in the same time 
a sufficiently current carrying capability together with the stability of the performances 
against electro-magnetic cycling loads. So, for the time being the critical current density 
J C (B,T) is computed according to the formula given in Ref. [13| . taking a nominal strain 
= -0.46%. 

The saturation regime model, and the evaluation of the condition for its applicability, 
additionally require the value of the empirical parameter -D e //. This quantity, in multistage 
CICC, is difficult to determine, although it places itself in the range between the filament 
diameter d e ff and the whole cable diameter [14J . In order to choose a reasonable value, we 
use as a first guess -D e // ~ 0.81 mm, the diameter of the actual SC strand, following [TS] . 



8 



The resistivity of the jacket is taken as p = 7.8 x 10~ 7 Qm, constant in the temperature 
range of interest, according to Ref. [TS]. The losses due to the eddy currents in the jacket 
are easily evaluated and result - as it could be expected - negligible, so that in the following 
we will ignore them. 

In order to set the temperature dependent quantities (J c for the AC losses, the densities 
and thermal capacities for the cooling part of the system), we take a constant and uniform 
T = 4.6 K. Obviously this is a crude approximation, but coherent with the scope of the this 
work, which is aimed to a preliminary assessment of the heat load in the CS superconducting 
cable. Refined simulations with thermo-hydraulic simulation codes such as Gandalf, taking 
into account the full time evolution and space distribution of the various physical quantities, 
will be performed in future works. 

The choice of the model for the AC losses computation (static, saturated regime, or 
crossover) should be done according to the saturation parameter value. A straightforward 
evaluation of the saturation condition (/3 > 1 — |i|) shows that the innermost turns go into 
the saturation regime at the beginning of the field ramp-down (see Fig. [2] for the innermost 
turn) . 

10 
1 



0.1 
0.01 

5 10 15 t(s) 20 

FIG. 2: Check of the saturation condition for the innermost turn: at t=0-4-0.04 s the saturated 
regime is expected. (Lines are guides for the eye) 

With the aim of yielding a preliminary evaluation of the AC losses, we choose to keep the 
used model as simple as possible, also in order to minimize the impact of a certain degree 
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of arbitrariness of the approximations implicit in complex models. The easiest solution 
would be to use the static model only, which would allow to have a worst case evaluation 
since, by neglecting the saturation occurrence, it overestimates the actual losses [HJ [33]. 
Nevertheless, deep in the saturation regime the overestimation can become very large: for 
example, with Eq. (II]) a P c « 40 kW at t = s on the first turn is easily obtained, instead 
of the ~ 3 kW given by the saturated limit expression and Q. In order to obtain more 
credible figures, we opt for the following choice. We use a piece-wise mixed model: each time 
that the local total power loss, evaluated in the static regime, is larger than the corresponding 
quantity within the saturation regime, the actual losses are taken as those given by the full 
saturation expressions. Obviously in this way the crossover region of pre-saturation remains 
still overestimated, but large overestimations are neverthless avoided. 

The plot of the total power losses in the selected against time is reported in Fig. [3j 

P(W) 
10 5 

10 4 

10 3 

10 2 

10 1 

-10 -5 5 10 15 20 25 t(s) 

FIG. 3: Total power loss, and contributions due to various mechanims, against time in the pancake 
subjected to the largest field. 

It can be seen that the charge and discharge phases contribute with a significant heat gen- 
eration, equal to 12.8 kJ and 16.1 kJ[T7], which is comparable with the one generated during 
the main phase. Since the latter is by itself large and challenging for the coil (see later), 
the charge/discharge phases (and the corresponding heat loads) will no further commented, 
assuming that some sort of softening (i.e. a reduction of the field sweep rate and an eventual 
additional phase with a steady field - which would allow the SC cable to recover from the 
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previously generated heating [IB]) would be possible in the definition of the scenario. 

We now focus on the main phase: the total energy dissipated through all the mechanisms 
in the whole — 20 s time interval is about 39.3 kJ, which is definitely a large value with 
respect to the cooling power presumably available. The latter can be roughly evaluated by 
computing the {mC) to t = 'Yl l i m i^i °f an the components of the cable (the NbsSn and Cu 
fractions of cable, the steel of the jacket and the helium flowing in the void fraction of the 
cable), where rrii and Cj are the total mass and thermal capacity (at constant pressure) of 
the i th component. Obviously the dominant contribution arises from the flowing helium, the 
physical properties of which (mass density and thermal capacity) are evaluated at 4.6 K, as 
already stated, and pressure 6 bar. The helium total mass is computed as the mass already 
present in the cable at t = s plus the additional contribution arising from a standard flow 
of 5 g/s lasting for 20 s (the duration of the main phase). One obtains (mC)tot=Q-7 kJ/K, 
definitely insufficient to remove all the heat generated by the AC losses. Obviously this direct 
computation does not take into account the actual time evolution and space distribution of 
the heat load and cooling power, which make the scenario even more challenging. Such a 
full analysis, which has to be made by thermo-hydraulic simulation codes - as e.g. Gandalf, 
is out of the scope of this preliminary work which, as already stated, aims to preliminarily 
estimate the heat loads and possibly indicate the general directions of the design choices 
and changes needed to satisfy the project requirements. 

In this perspective, it is worth noting that total energy Q t =39.3 kJ is the sum of Qh + 
Qd=l5.9 kJ and Q c =23.2 kJ (see also the energy dissipation distribution per turn reported 
in Fig. |]). 

This almost balanced splitting between coupling and hysteresis losses shows that, in order 
to significantly reduce the AC losses of the SC cable, both d e ff (impacting on hysteresis 
losses) and nr (impacting on coupling losses) have to be reduced. For example, by assuming 
an ideal extremely high-performance cable, with nr ~ 10 ms and d e ff = 10 /mi, one could 
reduce the total losses by a factor of 5, lowering them to 7.5 kJ. By taking advantage of the 
shortness of the coil length, one could choose a higher helium flow (8 g/s), which would yield 
10 kJ/K, so that the ratio Qt/(rnC)tot = 7.5/10 = 0.75 K is nearer to the needed operation 
conditions, but still unsatisfying. Indeed, the current sharing temperature T cs in the worst 
case, i.e. in the innermost turn at t=20 (so that one has the worst-case peak field 17.7 T 
and the maximum I t ), is ~ 5.3 K; if the above value of 0.75 K gives a uniform heating of the 
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FIG. 4: Total dissipated energy, due to the various mechanism, for each turn. 

SC cable, one would have 4.6 K + 0.75 K = 5.35 K > T cs , giving no temperature margin. 

Going back to the instantaneous quantities, in Fig. [3] it can be seen that coupling losses 
are dominant during the initial fast ramp down, whereas hysteresis losses are dominant in 
the flat top phase. The latter phase is the less challenging, although the dissipated power is 
still quite high, being on average ^194 W, whereas the cooling power is around 19 W/K. A 
possible improvement of the heat balance can be obtained with the same criteria as before: 
by taking d e ff = 10 /im (main intervention) and nr = 10 ms, losses drop to 40 W, whereas 
the improved cooling gives 31 W/K. In order to sustain the flat top phase, therefore, large 
modifications of the cable should be taken into account; perhaps, also an increase of the 
cable void fraction to increase cooling should be considered (compatibly with other design 
requirements). 

An even worst situation involves the initial ramp-down phase: an average 4.3 kW dissipa- 
tion occurs, with a peak of 140 kW in the initial 0.04 s, caused by the very high B involved. 
A plot showing the distribution of the various kinds of power losses over the 14 turns of the 
pancake separately is shown in Fig. |5j 

It can be seen that coupling losses are by far dominant. For comparison, the coupling 
losses computed within the static model only and the saturated regime model only are also 
reported. The wide maximum on P c around the 7 th turn is due to the switch from saturated 
regime to the non saturated one: the losses in the nearby turns are definitely overestimated 
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FIG. 5: Power loss at t=0 s per turn. 



because of the rough modelization here used. 

The roughness of the model obviously conspires against a sensible assessment of the actual 
margins: in these extreme conditions, in order to obtain trustworthy figures capable to guide 
the design process, both a full (and validated) model for the pre-saturated and saturated 
regimes is needed, as well as a reliable estimation of -D e //, which is the main scale factor 
governing both the entering in the saturated regime and the scale factor of the corresponding 
losses. Nevertheless, again an approximate consideration can be done: the total dissipated 
energy in the first 0.04 s is 5.6 kJ. Given the extremely short time interval considered, this 
quantity can be reasonably compared with the cooling power due to the helium already 
present inside the cable. Since the latter yields 4.7 kJ/K (in nominal conditions), the two 
figures would yield an increase of temperature ~ 1 K. By assuming the same ideal high- 
performance cable as before, hence by reducing by the usual 5 factor nr, d e ff and also D e ff, 
one obtains a smaller dissipated energy 140 kW/5 x 0.04 s=l kJ yielding an equivalent, more 
sustainable temperature increase of ~ 0.2 K. On the other hand, if one remains adherent 
to the presently available cable performances, this initial field sweep needs definitely to be 
slowed down: since the mission of the project dictated a given time derivative of the total 
B flux threading the plasma chamber, one could envisage a B reduction by simultaneously 
increase the radius of the CS. This solution, however, should be weighted against the already 
very high constraint due to the peak magnetic field on the conductor. 
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IV. CONCLUSIONS 



The AC losses heat load for the CS module for FAST, in the hypothesis of a supercon- 
ducting cable based design, has been evaluated within the very challenging H-mode reference 
scenario. This preliminary analysis shows that the high fields and high field rates regimes 
drive the SC cable to high dissipations conditions, both near and well within the saturated 
regimes. Therefore, it is mandatory to have full, reliable models, together with a proper 
characterization of candidate SC cables at saturation, in order to properly address these 
regimes and correctly estimate the corresponding expected heat load. Only in this way a 
proper evaluation of the real cable requirements can be done. 

Nevertheless, the present preliminary analysis allows to say that, in order to meet the 
stringent design requests, one has to consider much more performing strand (lower d e ff, but 
maintaining high Jc) and cable (lower nr and lower losses in the saturation limit). Even if 
some optimization of the cable design can be envisaged, most probably the field ramp rate 
should be reduced. More quantitative considerations will be provided in future works, in 
which full thermohydraulic simulations, based on the inputs of the present analysis, will be 
performed. 

In closing, we note that the difficulties of the present superconducting design are, at 
least partially, originated from the choice to rigidly adhere to the characteristics and volume 
constraints of the standard FAST design, which is based on resistive magnets. If it were 
possible to change completely the design in a superconducting cable perspective, maintaining 
in the same time the mission of the project, the requirements on the SC cable would probably 
become less stringent. 
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